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We study scalar cosmological perturbations in a braneworld model with a bulk Gauss-Bonnet 
term. For an anti-de Sitter bulk, the five-dimensional perturbation equations share the same form 
as in the Randall-Sundrum model, which allows us to obtain metric perturbations in terms of a 
master variable. We derive the boundary conditions for the master variable from the generalized 
junction conditions on the brane. We then investigate several limiting cases in which the junction 
equations are reduced to a feasible level. In the low energy limit, we confirm that the standard 
result of four-dimensional Einstein gravity is reproduced on large scales, whereas on small scales we 
find that the perturbation dynamics is described by the four-dimensional Brans-Dicke theory. In 
\^ • the high energy limit, all the non-local contributions drop off from the junction equations, leaving 

a closed system of equations on the brane. We show that, for inflation models driven by a scalar 

, field on the brane, the Sasaki-Mukhanov equation holds on the high energy brane in its original 

■ four- dimensional form. 
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I. INTRODUCTION 



Motivated by string theory, a new picture of our Universe has emerged, stating that our four-dimensional (4D) world 
is viewed as a "brane" embedded in a higher dimensional spacetime (the "bulk"). This "braneworld" picture offers us 
intriguing possibilities of testing theories with extra dimensions in future observations or experiments. Gravitational 
and cosmological consequences of braneworld models have been explored by a large number of references, and are 
\ reviewed, e.g., in 

One of the simplest realizations of braneworld is proposed by Randall and Sundrum (RS) 0,0], assuming that the 
\Q [ bulk involves five-dimensional (5D) Einstein gravity with a negative cosmological constant. The RS model can be 
naturally extended to include the Gauss-Bonnet (GB) term: 



Cgb ■= H 4 - m AB TZ AU + KabcdK a ^ u , (1) 

^ '" where 1Z, IZab, and Habcd denote the Ricci scalar, Ricci tensor, and Riemann tensor in five dimensions, respectively. 
This term arises in the low energy effective action of the heterotic string theory. The GB Lagrangian is the unique, 
ghost-free combination of quadratic curvature invariants leading to the field equations which contain derivatives of 
the metric of order no higher than the second Q- In the context of the brane model with the GB correction, 
linearized gravity in the GB braneworld has been studied in Refs. [H, @, 0], while nonlinear behavior of g ravity 
has been addressed [1, [§] using the geometrical projection approach of 10]. Cosmology on a GB brane [TTlfT^] is 
important as well, and one of the possible ways to test the braneworld idea is studying cosmological perturbations 
from inflation as they are linked directly to observations such as the cosmic microwave background. In this direction, 



Minamitsuji and Sasaki [13j have examined linearized effective gravity on a de Sitter (dS) brane, and Dufaux et al. 14 1 
investigated tensor and scalar perturbations generated from dS inflation in the GB braneworld (The authors of 14 ] 
have performed an exact analysis for the tensor perturbations, but they have neglected bulk effects for the scalar 
perturbations without any justification). In the present paper, we study scalar cosmological perturbations on a more 
general (flat) Friedmann-Robertson- Walker cosmological brane. 

Cosmological perturbations in braneworlds have been discussed in the vast literature, most of which focuses on the 
RS model and hence considers Einstein gravity [H, [H, E3, El, El, H3, [ML H2, HI, HI] • In this paper we basically follow 
the approach taken in the RS case, extending the previous results to include the effect of the GB term. 

This paper is organized as follows. In Sec. II we provide the field equations and the junction conditions in 
the presence of a bulk GB term and present the cosmological background solution. In Sec. Ill we consider scalar 
cosmological perturbations. First we derive the bulk metric perturbations in terms of a master variable, emphasizing 
that the 5D perturbation equations reduce to the same from as in the RS braneworld. Then we impose the junction 
conditions to give the boundary conditions for the metric perturbations. In Sec. IV we carefully investigate the 
limiting cases. Our conclusions are drawn in Sec. V. 
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II. GAUSS-BONNET BRANEWORLD 
A. Preliminaries 

We start with providing the basic equations that describe the GB braneworld. Our action is 



S = 2k 5 / d5x ^ - 2A + aL GB ] + J d 4 x 



4a 

2K+—Q + £ m -a 



(2) 



where A is the cosmological constant in the bulk, C m is the matter Lagrangian on the brane, and a is the brane 
tension. The GB Lagrangian £gb was already defined in Eq. ([T]) and the coupling constant a has dimension of 
(length) 2 . The surface term is given by 2K + (4a/3)Q, where K is the trace of the extrinsic curvature of the 
brane and Q := Q £ with Q" defined below in Eq. ([B]). 

The 5D field equations following from the above action are 

Gab — ~2^-ab — —h.9AB, (3) 
where Gab '■= T^-ab ~ TZ-9 ab/2 is the Einstein tensor and TIab is the GB tensor defined by 

Hab ■= C GB g A B - 4 (TZTZab - 2H AC Tl c B - 2H A cbdH CD + Hacd E iiE de ) . (4) 
Assuming a Zi symmetry across the brane, the junction conditions at the brane are given by [25j ] 

V - KS; = ~ {T; - aS;) - 2a (q» - Iqsj) , (5) 

where T^ v is the matter energy-momentum tensor and 

Q; := 2KK«K£ 2K«K c ?Kf + {K^Kf - K 2 ) K» 



2KR; + RK; - 2K£R^ a - 2R«K a » - 2R a v K«, (6) 



with R^va/3, R^v and R being the Riemann tensor, Ricci tensor and Ricci scalar with respect to the 4D induced 
metric. The main difference in the junction conditions from those in Einstein gravity is that they include intrinsic 
curvature terms as well as external ones. As we will see later, it brings significant modifications to the behavior of 
cosmological perturbations. 

Using the Codacci equation, we can show that the conservation law holds on the brane @ : 

V U T^ = 0. (7) 

The field equations admit an anti-de Sitter (AdS) bulk with the curvature radius I (—: The 5D cosmological 

constant and /x are related by 

A = -6fi 2 (1 - 2 af i 2 ) . (8) 

It is useful to define a dimensionless parameter (3 := Aafj 2 . In this paper, we assume the parameter range < /3 < 1. 
The upper limit here is sufficient for ensuring A < 0, but as we will see later the linearized field equations are associated 
with the overall factor (1 — and therefore we assume the tighter limit. 

B. Cosmological background solution 

We present a cosmological background solution which has a flat 3D geometry in the GB braneworld [ll| . We write 
the metric in the Gaussian normal coordinates as 

g%dx A dx B = -n 2 (t, y)dt 2 + a 2 (t, y)5 i3 dx i dx 3 + dy 2 . (9) 

We may set n(t, 0) = 1, so that t is the proper time on the brane at y = y\, = and a&(t) := a(t, 0) is the scale factor. 
The 5D field equations for this metric are given in Appendix [a1 If Aab ^ 1, where 

>:=(£V-i(iY, (10, 
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where a prime (overdot) denotes a derivative with respect to y (t). Eq. (|A3|) requires (a/n)' = 0, and hence 

n{t,y) = ——. 11 

Substituting this into Eq. (|A1|) . we obtain [a 4 6 — 2aa 4 6 2 ]' = —(A/6) (a 4 )' , which can be integrated immediately to 
give 

b-2ab 2 = p 2 (l- ^) +4- ( 12 ) 



The integration constant C corresponds to the non-zero components of the Weyl tensor in the bulk. In this paper, we 
restrict our analyses to the AdS background and assume C — 0. For vanishing C we obtain 

b = p 2 . (13) 

We do not consider another possibility, b = (2/f3— l)p 2 , because this does not provide a well-behaved a — ■> limit. 
Evaluating Eq. (fl3| at the brane, we obtain 



= -vV + # 2 , (14) 
where _ff := db/at, is the Hubble parameter. Then Eq. (|11| implies 

< = -^WT7- f (15) 

Eqs. (|14[) and (|15[) will be used in the next section when discussing cosmological perturbations on the brane. Substi- 
tuting b — p 2 into Eq. (|A1[) , we obtain a"/a = p 2 . From this and Eq. (TTJJ we find 



a(t,y) = a b (t) 



/ # 2 
cosh(py) - J 1 + — sinh(//y) 



(16) 



Although the 5D field equations include the GB term, the metric functions n(t, y) and a(t, y) have the same form 
as in the cosmological solution in the RS braneworld based on the Einstein-Hilbert action [26|. What is manifestly 
different is the Friedmann equation that relates the Hubble expansion rate H and the energy-momentum components 
on the brane. The Friedmann equation derived from the generalized junction conditions at the brane is [l2| 

2^H 2 + p 2 [i - 13 + W^P) = « 2 (/> + (17) 

The critical brane tension, which allows for a Minkowski brane, is obtained by setting H — * as p — > 0: 

k 2 ct = 2^(3-/3). (18) 

There are three regimes for the dynamical history of the GB brane universe, two of which are basically the same as 
those found in the context of the RS braneworld. When H 2 <C /i 2 /P[= (4a) -1 ], we recover the RS-type Friedmann 
equation, 

, 8^G ( _ , p 2 



— (19) 
where we defined the 4D gravitational constant as 

S. G :=^. (20) 

Thus, we can see that at low energies, H 2 -C p 2 , we have the standard 4D Friedmann equation, H 2 cx p, while at 
high energies, p 2 <C i? 2 (<C p 2 //?), we have H 2 cx p 2 . At very high energies, H 2 ^> p 2 /(3, the effect of the GB term 
becomes prominent. In this regime, we find 



4/? 
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III. COSMOLOGICAL PERTURBATIONS 



A. Perturbations in a maximally symmetric bulk 



Now let us consider linear perturbations about the cosmological brane background discussed in the previous section. 
Since we are considering the maximally symmetric bulk spacetime, the background Riemann tensor can be expressed 

as R^bcd = — A* 2 9ac9bd — 9ad9bc • Using this fact, we find that the perturbed GB tensor has a following nice 
property: 

Thus, the linearized field equations are simply given by 

(1 - p)5Qf = 0, 



(22) 



(23) 



which, aside from the factor (1 — /3), give the same perturbation equations as in Einstein gravity. This allows us to 
make full use of the previously known results on cosmological perturbations in the RS model. Note that Eq. (|22| is 
a direct consequence of maximal symmetry of the background. 
We write the perturbed metric in an arbitrary gauge as 



[9% 



2E tij ] dx l dx 3 



5g AB jdx A dx B = -n 2 (l + 2A)dV + 2a 2 B^dtdx l + a 2 [(1 

U JV )dy 2 . (24) 

The gauge dependence of the metric perturbations is summarized in Appendix [Bj The 5D perturbation equations 
will be solved most easily in the so-called 5D longitudinal gauge [l6j . which is defined by 



a = -B 

7y By 



E = 0, 
f E' = 0. 



(25) 
(26) 



(Hereafter variables with tilde will denote the metric perturbations in the 5D longitudinal gauge.) We use a master 
variable, Q, which was originally introduced by Mukohyama [IH in the Einstein gravity case. The perturbed 5D field 
equations are solved if the metric perturbations are written in terms of this master variable: 



A = - — 



4> 



1 

6 a 
1 

na 
1 

6a 

1 

6a 



2n" - -ft' 

n 



2 n 



4 1 <> 



n" 



— n 

n 
n 



n ■ 

-n 

n 



-n 



n" 



-o' - 2n 2 n - 



where Q is a solution of the master equation 



n" 



, n 

n - | - 

n 



3- ) n 

a 



~2 A 



n = 0, 



(27) 
(28) 
(29) 
(30) 

(31) 



with A := S^didj. The boundary conditions for will be derived in the next section. 

The master equation (f3"Tj) does not have a separable form except for the special case of a dS brane background. (The 
separable dS braneworld is discussed in Appendix[Cl) As has been worked out in the context of the RS braneworld [24| 
(see also [2l|, [22|, [23j]), one must in general resort to numerical calculations to solve Eq. (f3"Tj) . 



B. Junction conditions 



Perturbed junction conditions provide the relation between the gravitational and matter perturbations on the brane, 
which leads to the boundary conditions for the master variable. The junction conditions are most easily derived in 
the brane- Gaussian normal (GN) gauge [IB], which is defined by 

(9% + S ~9ab) dx A dx B = -n 2 (l + 2 A) dt 2 + 2a 2 B^didx l + a 2 [(l - 2$) + 2E tij ] dx l dx j + df , (32) 
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and 



y b = o. 



(33) 



Here we denote by a bar the perturbations in the brane-GN gauge. Starting from the 5D longitudinal gauge in which 
the location of the brane is perturbed and is given by y b = the brane-GN gauge is realized by a coordinate 

transformation x A = x A + Sx A such that 







1 



B y - Sx' -Sy 



A v + nSt' — —Sy, 
n 





= Ayy - 5y', 



and 



= £ + Sy b . 



(34) 



(35) 



There is a residual gauge-freedom in the time coordinate, and in the following discussion we fix 5t b = 0. 

We decompose the spatial component of the perturbed extrinsic curvature into its trace SKt and traceless part 
SKtl as 



SK T S, J + 



SK- 



TL- 



In the brane-GN gauge, the extrinsic curvature is simply calculated as 

SK ° = A', 
SK,° = 



1 2 1- -i\ 
-a b (a y -a ) t , 



SK T 
SKtl 



-j'+^ACTy, 



(36) 

(37) 
(38) 

(39) 
(40) 



Hereafter in this section all the perturbation variables are evaluated at the brane. A straightforward computation 
gives 



SQo° = -G^ 2 SK ° - 12l?8K t + 2 



a b n b 



SG ° - 6-±6G T , 
a b 



SQ t ° = -2^5^ + 2^50°, 
a b 



SQt 
SQtl 



2 fa' h 



4fi 2 SK ° - U^ 2 5Kt -rA- + 1-)5G ° - 4—SGt, 



3 \a b 

-2^SK TL + 4 ( _ 



m, 



a b 



a b 



(41) 
(42) 
(43) 
(44) 



where the trace and traceless part of SQ/ and the trace of SG/ are defined similarly to SKt and SKtl ■ The perturbed 
4D Einstein tensor is given in terms of the induced metric perturbations as 



SG ° = 6hU + HA) - 4 A *> 
SG t ° = -2($ + HA)., 



SG t = 2 



ip + 3Hip + HA + (3H 2 + 2H) A 



2 1 
3^ 



2 A(*-$). 



(45) 
(46) 
(47) 



and the metric potentials are defined by 



A- S (afc), 
i\) + a 2 b Ha. 



(48) 
(49) 
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The perturbations of the energy-momentum tensor are given by 



oT ° = -5p, 
ST, = Sq ti , 



,$T/ J = SpS/ + ( 8$ - ~5f A ) Sir. 



(50) 
(51) 

(52) 



From the junction conditions {5| we obtain 



K 2 S P = -6(1 -P)8Kt + ^^SG °, 

K 2 % = -2(1-/3)^°-^^°, 

M «6 



K 2 (5p = 2(1 - /3) (^ ° + 2<5At) + 



k 2 Stt = -2{l-p)8K TL -%\ 





1 








3 




rib J 














a b 





SG ° - ^SG T 
ab 



(53) 
(54) 
(55) 
(56) 



It is worth noting here how the perturbed 4D Einstein tensor appears in the junction conditions. Eqs. (|53p and ([54 
clearly have the suggestive form of 



k 2 5T^ = (extrinsic curvature) - ^f^<5G/. 

\X (X}) 



(57) 



The other two equations (|55|) and (|56|) have the form slightly different from (|57l) . but when a' b /ab — n' b /rib oc H = 
they reduce to ((57)) . 

In order to write the junction equations in terms of the master variable O, we now go back to the 5D longitudinal 
gauge by substituting Eqs. (|34|) and (|35|) into the the gauge transformation (|B2|) . The extrinsic curvature is expressed 
in the 5D longitudinal gauge as [16J 



SK n 



A + A>, — -A vv + £ 

"6 



i 

2a b 



3a? 



A ( fi' - — fi 





-<) 


2< 


-<) 




rib) 


V a b 


rib) 



i 

2ab 



2H 



< 


<\ 


K 


a b 


rib) 


rib 



rib 



5JC T = 



<5^TL = 



1 

2afc 



f2' 



rib 



-& + HA y -^A yy + H{t-HZ)-~AZ 
ab 



1 1 



2a b 







6a? V a 6 



3 a 2 



1 1 

3 a 



2a b \a b n b 



n 



(58) 

(59) 

(60) 
(61) 



Similarly, in terms of the 5D longitudinal gauge perturbations the 4D Einstein tensor and the metric potentials are 
expressed, respectively, as 



SG V = 6H [rp + HA -— Atp - 6 



. a,. 



at 



ab 



a b 



b a b 

/i 2 + H 2 + H)il 



AC 



i, 



a 



1 1 

^n' - U 2 + h 2 ) n 

a b 



Za 5 b 



A 2 n 



7 



-6- 



06 



HM-Ht)-~A£ 



SG° = -2 



SGn 



jj + HA-^ (i- HA 
a b V 1 





~<h> ( 






yah 


a b \ 



H + H)fl 



3al 



A n - HQ. 



2-M£-if£ 

a b 



(62) 



(63) 



i/> + 3Hi> + HA + 3iJ 2 + 2H) A 



a b \ a b n b 



2 1 
1 1 



A 



1 

a b 



^n' - ( M 2 + h 2 ) n 





K 


I- ^ fffi' - 


l r 
f — 




\ a 6 




a b L 



1 

a b 



2^H 

a b 



^H 

n b 



a b 



n b J 3at 



2^ ( % +2 ^W 

2H (^ 2 + i? 2 
1 



Sl'-\ I n' 

a b 



-2^| 
a b 



(<- 




\a b 


n b J 



3og 



ffA 



6.H\ff + i? 
1 



3 at 



2ff 2 [ " b 



9a£ 



a^ 



2^ U, 
ab J 



and 



4> 



a;, 



The matter perturbations are subject only to a temporal gauge transformation, and since we fix 6t b = 
invariant when going from the brane-GN gauge to the 5D longitudinal gauge. 
Looking at Eqs. ([55 ]) -([50 1) and (|52 l) -([53 |) carefully, we find 1 



1 (a' 



a b n b 



6G ° 



^6G ° 
a b 

^SG° 
a b 

a',. 



1 



} -5Gi 



-6 U 2 + H 2 ) SK T + ^/p 2 + H 2 -^A 2 n, 
2 ( M 2 + H 2 ) 5K» + vV + tf'^A (p - HQ) 
2 (fi 2 + H 2 ) (SK ° + 2SK T ) 

H \ 1 



+AH6K T + vV + H 2 1 - 



H 2 + H 2 I 9al 



A 2 n. 



Using these equations, the junction conditions can be rewritten as 

M 2 



n z Sp = -6(1+/? + 213'— ) SK T 



2/3 



3/x V [i 2 a\ 



K 2 Sp 



M 2 
ft 

T2 



n 2 8q^ = -2(1 + /3 + 2/3^- )SK i ° - ^Wl + ^-^Afn-ffS) 



3[i 



H 2 1 

Ai 2 al 



AHSKt + vV + H 2 1 



if 



H 2 ) 9al 



A 2 n 



(64) 

(65) 

(66) 
they are 

(67) 
(68) 

(69) 

(70) 
(71) 
• (72) 



The traceless part of the junction equation (p)6"|) with Eqs. (|6"Tj) . (|65[) . and (j6"6")l yields 



l + /3 + 2-^ ( // " 



e = £ a i 5n+ JLL\<n-i2^-^)Hn + 2^[^-^)si' 



2fi 2 a b I a b 



a b n b 



a b \a b n b 



1 This is a consequence of a perturbation of the contracted Gauss equation. 
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At" — - — 

a b n b 



2°± _ Hi 
a b n b 



(73) 



In the RS braneworld (/3 = 0) , the brane bending in the 5D longitudinal gauge vanishes in the absence of the matter 
anisotropic stress, St: = 0. In the GB braneworld, however, Eq. (|T3|) clearly shows that the curvature tensors in the 
junction equations act as an anisotropic stress source, and hence the brane bending £ should be taken into account 
in general even if 5ir = 0. The situation here is in a sense similar to that of the Dvali-Gabadadze-Porrati (DGP) 
braneworld 27]: a induced gravity term on a DGP brane mimics anisotropic stress (see, e.g., [28j|). 

Using Eq. ([75]) (with Sir — 0), we can eliminate £ in Eqs. I|58 p — (|60 p . We then specify the relation between Sp and 
8p, for example, in the form of the equation of state Sp — wSp. This procedure leads to the boundary conditions for 
the master variable fl, which are in general quite complicated. 

A comment is now in order. Using the contracted Gauss equation [Eqs. f6"T |) -([59" |) ]. we can express the junction 
conditions in terms of the perturbed Einstein tensor rather than the extrinsic curvature. However, the junction 
conditions cannot be written solely in terms of the Einstein tensor on the brane; n itself will be inevitably involved. 
Such extra terms correspond to the 5D Weyl tensor, and for this reason closed equations on the brane are not available 
in general. 

IV. LIMITING CASES 

In the previous section we have presented our general formalism for scalar perturbations in the GB braneworld. We 
now investigate limiting cases in which the junction conditions can be simplified to a great extent, in order to catch 
on to the specific new features brought by the GB term. 



A. Low energy limit 



First let us consider the low energy limit maxjiJ 2 , \H\} <C p 2 (< [i 2 /f3)- In this limit, we may approximate 
a' b /a b ~ n' b /n b ~ —p and p (a' b /a b — n' b /n b ) ~ H. Then the junction equations (|70p - (|72p read 



n 2 Sp 
n 2 Sq 
K 2 Sp 



-6(1 + /?) 
-2(1 + (3) 
2(1 +0) 



1 

2a& 
1 

2a, 



2/3 1 A 2 



3p a b b 



2/3 1 
3p a?' 



AO- HCl 



2a,b at 



1 HF+£ + 2Hi-^A£ 



2a b 



3a 



2/3 1 , 2 



9p al 



— — | -<) -- //<> 

2p a b 



1 



An 



where 



^(t,x) : = [n' + pn] 



b ■ 



(74) 
(75) 
(76) 
(77) 

(78) 



and we set Sir = 0. We have dropped terms such as H 2 ^ and HI; in the above because Eq. ([77]) implies H 2 ^,H^ <C 
{p 2 / (3)£, ~ HJ- ja b . Defining a new variable 



V(t,x) := (1 + /3) {T + 2oi0 + 



2JS_ 

p, 



Hn + — An 

3a 2 



= (1 + /?)^- 
we obtain the following simple set of equations: 

n 2 a b 5p = -3HV 



n - 3Hn - —An 

a 2 



jAV, 



K 2 a b 5p = V + 2HV + HV. 



(79) 



(80) 

(81) 
(82) 
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1. Perturbations larger than the bulk curvature radius 



Now we assume (as in 19]) 



fj,n', n 2 n » n, ho,. 



Then we have 



JL 

2ab 

JL 

2a b 



1 2 + (3 
6a? 1 + f3 



An, 



T + ^-±^An, 

6^ 1 + /3 



and 



V « (1 + /3).7 r + -^Afi. 

On scales much larger than the bulk curvature scale, the gradient terms in Eqs. 
\\iT\ ~ \[i 2 n\ ^ AQ/af . Consequently, we obtain 



1 



r A* = AttG Se, 



(83) 

(84) 
(85) 

(86) 

can also be neglected as 

(87) 
(88) 



and 



$ + (4 + 3c 2 ) iJ$ + [2JT + 3 (l + c 2 s ) H 2 ] $ - c 2 ^ A$ = AnG (Sp - c 2 s Sp) 
L J a b 



(89) 



where 5e := 5p — 3H5q is the comoving density perturbation and c s is the sound velocity. Thus the standard 4D result 
is reproduced. Since Eqs. ([56]) and (j6"Tj) with vanishing matter anisotropic stress imply 



$ - * 



1-/3 




(90) 



there is no brane bending in this case. 



2. Small scale perturbations 

Let us turn to scales much smaller than the typical GB scale, |Af2/a 2 | ^> |/Lt 2 r2//3|. For the moment we also assume 
the quasi-staticity of ft [Eq. (|53")) ] 2 . In this case we may ignore F terms relative to gradient terms in Eqs. (f8"4"|) - ([56]) . 
so that 

.La* = 4r (2+2) g &, (9 i) 
-* A ,_«,(I±*) Gfc (92) 

These are the cosmological extension of the result in Q . In this case we have $ — W ^= and hence the brane bending 
plays an important role. Since we are considering length scales smaller than P 1 ^ 2 /^, we cannot take a smooth limit 

The above result implies that the perturbation dynamics on small scales is described by a scalar-tensor type 
theory We can show that this is in fact the case without invoking the approximation (|83p . Here we only require 



2 Note that this approximation is valid only for dust matter because on small scales we have fl ~ (c 2 fc 2 /a^)S7 in general. 
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£ ~ (0/[i)A{l/a%. This in particular means that f2 can be as large as O(A0,/ a 2 ). Then our approximation allows one 
to neglect fi terms (which can be expressed solely in terms of T at low energies) relative to the brane bending £ in 
the perturbed extrinsic curvature (f58"|) - ([60"l) . As a result, Eqs. (|53|) -([55 |) can be written as 



SG Q ° 



«V (1-0 



8G1 



k 2 h 



1-0 



1-/3 




3^-Ia^ 

af 



where we used A£/a 2 3> iJ 2 £. Neglecting terms in Eqs. ([74]) and ([76]) . we also have 



6(1 + 0) 



{-Sp + 3S P ) 



6(1 + 0) 



ST 



The above four equations and Eq. (|90[) are equivalent to 



= ^-st; + — (v M v 

" 2(^o A ¥>o 
1 



V A V A <5/) ty, 



(93) 
(94) 
(95) 

(96) 
(97) 



6 + 4w 



ST, 



with the identifications 



1 



</?o 



1-/3 
/3 



3/3 
1^0- 



(99) 



This is nothing but the linearized Brans-Dicke (BD) theory with terms of 0(H 2 5ip) neglected. The result here is in 
agreement with the previous results for a Minkowski brane 0] and for a dS brane in the low energy limit (T3 |. 

In Ref. @ it is argued by looking at Newton's potential that in the GB braneworld one can take the scale £ (— p^ 1 ) 
to be of geophysical size 3 (say ~ 1 km - 100 kms) with not too different from unity (say > 0.85). One would take 
£ to be much larger and at the same time fine-tune to be extremely close to unity, so that the BD-type theory 
would pass astronomical tests in the Solar System. However, in order not to spoil the standard cosmology picture 
after Big-Bang nucleosynthesis, I must be smaller than the Hubble horizon size at that time. Therefore, £ is required 
to be < O(10 13 cm) ~ O(AU) and so must be below cosmologically interesting length scales. 



B. High energy limit 



Now let us take the high energy limit, H 2 3> p 2 / 0. In this limit, the junction equations in the form of (|53 p ~ (f5"6")l 
are more convenient. Eq. (|67p implies 



p Z CLb 



» \SK T \. 



(100) 



We have assumed that the first and second terms in the right hand side of Eq. (|6"T|) do not cancel each other. This 
shows that the right hand side of the junction equation (|53[) is dominated by the perturbed 4D Einstein tensor. The 
same is true for the other two equations (|54p and (|55|) . Thus, in the very high energy regime we have 



SG» = 



SG,° 



SGt — 



20H 



Sp, 



20H 
20H 



Sq,i, 



(s P - e fsp), 



(101) 
(102) 
(103) 



3 This is in contrast to the RS model in which the bulk curvature radius is constrained to be I < 0.1 mm. 
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where we defined 

*n'=-^2' ( 104 ) 
Similarly, it is easy to show that the right hand side of Eq. |56[) is dominated by the metric potentials. Thus we obtain 

2 2 

(1 - e H )* - * = yjgal&it. (105) 

The set of equations (|L0ip — (|103|) and (| 1 05[) governs the perturbation dynamics at very high energies, H 2 ^> n 2 / (3. 
(In this regime, we cannot take a smooth limit /3 — > 0.) We should emphasize that all the non-local terms drop from 
the junction equations and hence the system is closed on the brane. Consequently, we know about the evolution of 
perturbations without solving the bulk. We remark here that these effective equations are shown to be consistent 
with the Bianchi identity by using the background equation (|2ip . 

For a dS brane background (ch = 0), which is the situation studied in Ref. the above equations reduce to the 
perturbed Einstein equations, 

SG; = 8ttG cS (106) 
where the effective gravitational constant is given by 

to** := %■ (107) 

Thus, recovery of Einstein gravity on a dS brane in the high energy limit (l3j is confirmed. For a general background 
with en =/= 0, the perturbation equations differ from the Einstein equations and the effective gravitational coupling is 
time-dependent . 

Now let us investigate the case with H =/= in more detail. Consider braneworld inflation driven by a single scalar 
field <j) which is confined on the brane. For this background we have p = <f) 2 /2 + V(<f>) and p = <p 2 /2 - V(<f>), where 
V{4>) is the potential of the inflaton. For perturbations generated by fluctuations of the scalar field, it is quite e asy 
to describe the evolution of perturbations in the high energy limit by introducing the Sasaki-Mukhanov variable [291 ] 
and invoking the energy conservation equation. The perturbations of the energy-momentum components are given 

by dp — 4> (s'(f> — <f>Abj + (dV/d(f>) 8(f), Sq — —$8<j>, and bp = <j> (5(f) — <j>Abj — (dV/dcj>) 5(f). The equation of motion for 

the scalar field perturbation 5(f) follows from the energy conservation equation, 5 (V V T^ U ) — 0. In terms of a scalar 
field perturbation in the spatially flat gauge, 



-g^b, (108) 



we obtain the wave equation 



*» + art*, + (f + 0)*-^4 (IS) «** = »■ <™) 

Here perturbations were Fourier decomposed as usual, with k being comoving wave number. In deriving the wave 
equation we used the field equations (|101[) - (|103p and (|105|) . and hence the last term looks different from the corre- 
sponding equation in standard 4D cosmology. However, introducing new variables 



v := ab5(f>,p, z := (HO) 
H 

Eq. (I109P can be rewritten in a familiar form 

v"+(k 2 -j\v = 0, (111) 

where a prime denotes a derivative with respect to conformal time r\ := J a~ 1 dt. This exactly coincides with the 
Sasaki-Mukhanov equation derived in the standard 4D context [29l |. Since the comoving curvature perturbation 
1Z C := v/z is conserved on super-horizon scales irrespective of the gravitational field equations [30| . it is natural that 
the conservation equation can be recast into the form of v" ~ (z" / z)v as k 2 — > 0. In the present case, contrary to 
the RS model, the brane and bulk perturbations are decoupled on small scales. Therefore, our result might not be so 
surprising. Even so, we believe it interesting enough to emphasize. 
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V. CONCLUSIONS 

In the present paper we have presented a formulation for scalar-type cosmological perturbations in a braneworld 
model with a bulk Gauss-Bonnet (GB) term. As a background solution, we have considered a 5D anti-de Sitter (AdS) 
bulk with the curvature radius (JT 1 , bounded by a flat Friedmann- Robertson- Walker cosmological brane. We have 
assumed that the AdS curvature radius is larger than the typical length scale defined by the GB coupling a (more 
precisely, (3 :— Aa/j 2 < 1). We have also assumed an arbitrary expansion rate H of the brane universe. The bulk 
GB term does not change the 5D perturbation equations from those in the Randall-Sundrum (RS) model, and hence 
we have adopted the approach using Mukohyama's master variable which was first introduced in the Einstein gravity 
case. As for the boundary conditions at the brane, the generalized junction equations in the presence of the GB term 
bring several new terms, which can in principle cause dramatic changes of the behavior of perturbations in comparison 
with the RS braneworld. 

In order to clarify the effects of the GB term, we have investigated the limiting cases where the boundary conditions 
are simplified to some extent. We have shown that in the low energy limit, H 2 <C /i 2 , the RS model is reproduced 
on large scales, a^X 3> M > where A is the comoving wavelength of perturbations. Namely, gravity on the brane is 
basically described by general relativity in this regime. On small scales, however, the result is quite different from 
that in the RS model. We have shown that the behavior of perturbations is effectively governed by the linearized 
Brans-Dicke theory for a^A <C /J 1 / 2 /^. 

At very high energies, H 2 3> /x 2 //3, the presence of the GB term leads to the most significant changes. We have 
found that in this high energy limit the evolution of perturbations on the brane can be determined without reference 
to the bulk perturbations. This is because the perturbed extrinsic curvature terms in the junction equations are 
suppressed compared with the novel terms arising due to the bulk GB correction, the latter being expressed solely in 
terms of the local quantities on the brane. 

Finally, we shall comment on scalar perturbations generated by the fluctuations of an inflaton on the brane. In order 
to determine the amplitude of scalar perturbations, one has to quantize the perturbations. In the RS braneworld, the 
perturbations on the brane are strongly coupled to the bulk metric perturbations on small scales [20l Izij , and hence 
the quan tization of the coupled brane-bulk system is required. This is the outstanding challenge in the RS case (see, 
e.g., [3l|). In the high energy regime of the GB braneworld, however, the perturbations on the brane and in the bulk 
are decoupled irrespective of their wavelengths, which enables us to quantize the system straightforwardly. Thus our 
result validates the approximation of 14] ignoring bulk effects only in the high energy limit. As the energy scale of 
the brane becomes lower (H 2 < fi 2 /(3), the interplay between brane and bulk perturbations becomes strong for short 
wavelength modes, and one will face the same problem as in the RS braneworld. 
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APPENDIX A: BACKGROUND FIELD EQUATIONS 



The 5D field equations are given by 
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Go" ?W„° = 3 
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where 



and an overdot (prime) denotes a derivative with respect to t (y). 



APPENDIX B: GAUGE TRANSFORMATIONS 

Under a scalar gauge transformation, 

t -> t + St, 
x i -> x l + d^x, 



the metric perturbations transform as 



A -» A- 5t--5t-—Sy, 
n n 

a c a ' c 

ip — > v + H — oy, 

a a 
n 2 

B -> B+-5-<ft-&E, 

or 

B y —> B y - 6x' 5-<5y, 

ar 

A y — > A y + n#t' <5y, 

j4 y!/ — » — 5y . 

It is useful to introduce the spatially gauge-invariant combinations 

a := -B + E, 
a y := -By+E, 
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which transform as 



n 2 



g — ► ct — —5t, (B5) 
<j y — * cry + \5y. (B6) 

APPENDIX C: INFLATON PERTURBATIONS ON A DE SITTER BRANE 

In this appendix, we consider a inflaton field <j> whose potential V is very flat: dV/d<fi ~ 0. For such a potential, we 
may assume a de Sitter (dS) brane background. Although the general discussion of gravity on a dS brane with the 
GB correction was already given in [l3| , for completeness we shall revisit the issue here using the formalism developed 
in the main text. 

Our approximation is as follows [2(| (see also 28]). First, we take the slow-roll limit and work in the dS brane 
background. We also assume that we can neglect the brane metric perturbation contributions to the matter pertur- 
bations, which is a valid approximation in the slow-roll limit of the standard 4D calculation in the longitudinal gauge. 
Of course, such a simplified description of the inflationary universe may be a toy model, but this is a price to pay for 
a feasible problem; the master equation (|3ip is separable for a maximally symmetric brane and hence we are able to 
obtain a bulk solution for the master variable $1 analytically. 

In the dS brane background, we have 

n(t,y)=n(y) = cosh(/xy) - 7sinh(^ty), (CI) 
a{t,y) = a b (t)n(y), (C2) 



where for notational convenience we defined 7 := yl-j- H 2 //j 2 . There is a Cauchy horizon at y = yu = H 1 coth 7. 
Now it is clear that the master equation (f5T|) is separable for this background: 



n" - 2— n ; + v 2 n - \ 



0. (C3) 



We write the solution to this equation in the form of 

fi(t,y,x)= / d 3 kn k (t, y y k *, n k (t,y) = I dm <p m (t)xm(y). (C4) 



From now on we will work in the Fourier space and suppress the subscript k. The mode functions satisfy 

ip m - SHtp m + (m 2 + ^ ] ip m = 0, (C5) 

- 2- X ' m + ( — + M J Xm - 0. (C6) 

In terms of the conformal time rj := — l/(abH), the general solution in the time direction is given in the form of a 
linear combination of Bessel functions of order v: 

Vm = (~fc^r 3/2 [ Cl (m) M-krj) + c 2 (m)Y u (-kr,)} , v 2 = - - ^J, (C7) 

4 H z 

where ci(m) and C2(m) are constants. The general solution in the extra direction is obtained in the form of 

Xm = n(y) [c 3 (m)P u _ 1/2 (coth fi(y h - y)) + c i {m)Q v _ 1/2 (coih h i(y h - y))] , (C8) 

where P a and Q a are and Legendre functions of the first and second kind, of order a respectively, and 03(771) and 
04(771) are constants. 
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The junction conditions are reduced to 



K *5 P = -6(1 + /?) [^.HF + ^F + Hi-^ + ^ 
,-r<u, = -2(1 + ( ±-F + t-HZ 



3/x a a b 



2/3 7 k 2 



1 ^. 1 „^ i- _„i _„o. 2fc 2 



ft &p = 2{1 + P)[ —F + —HF + t + 2Ht- 3H 2 C + 



2/i a;, 



2a fc a & 

1.2 

3af 



2f3j k 4 
9/i a 



n. 



5 uo ' 
6 



where 



and 



Dehning 



F(t) := fi' + /ryfi, 

(3 := -/i + 2/3 7 2 
= (3 + 2(3H 2 /n 2 . 



V ds (t) := (l+ / S)(F + 2a 6 + ^ (fffi - H 2 !! - Ajfl) 



07 



Q - 3Hil + 2H 



k 2 



Eqs. (|C9|) - (|C11[) are rewritten in the following simple form: 

VdS = k 2 
V ds + 2HV ds = n 2 a b 



dV 



dV 



(C9) 
(CIO) 
(Cll) 
(C12) 

(C13) 
(C14) 



(C15) 

(C16) 
(C17) 
(C18) 



where the energy-momentum components are now given by the perturbed scalar field. Combining Eqs. (|C16[) - (jC18[) 
and using the background equation cj> + ZH<j) + dV/dcj) — 0, we obtain 



PdS- ( H+t\V dS + ^V dS = 0. 



(C19) 



Keeping in mind the slow-roll condition \4>/4>\ ^ the solution to this equation is given by 

cos(-fc??) sin(-fc?/) 

/ (IS = t-i ; r O2 : , 

—krj —krj 

where C\ and C2 are constants. This solution can be expressed as a sum of Bessel functions as [2C 



(C20) 



T dS = C^f^i-l) 1 (21 + 1) (-k v )- 3 / 2 J 2l+1/2 (~k v ) 



1=0 



+c 2 x/2^(-i) ; (21 + -) (-M" 3/2 ^ +3 / 2 (-M- 



(C21) 



This indicates that the boundary condition is satisfied by an infinite sum of discrete mode solutions (— krf) 3 / 2 J v {—krj) 
with v 2 = 9/4 — m 2 /H 2 , where the mass spectrum is given by 



m 



-2{2l - + 1)H 2 for Ci, 
-21(21 + 3)H 2 for C 2 . 



(C22) 
(C23) 
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Let us construct the bulk solution consistent with the boundary condition. Due to the boundary condition induced 
by the scalar field on the brane, we can choose only the normalizable modes Q a in the extra direction. Then the the 
bulk solution is given by 

nfa y ) = Cifiifa, y) + c 2 n 2 (v, y), (C24) 

with 

n ( \ FT^t W /oi i X \ n(y)Q 2l (coth^(y h -y)) 3/2 

^ V 2 J H [(1- (3)Ql l+1 (j)~ {H/fi)PjQ 2 2l+1 (j)\ 

where Q^ a are associated Legendre functions of the second kind. 

Since the mode solution in the time direction behaves as (— ki])~ 3 / 2 J„{— krf) ~ a b ly+3 / 2 n super-horizon scales, the 
dominant mode is the one with m 2 = 2H 2 (the C\ mode with 1 = 0). Thus, on super- horizon scales we have 

n « C iab (t)tx(y h - y)n(y), (C27) 

where unimportant factors are absorbed into the redefinition of C\. Using this fact and Eq. (|30[) . we obtain 

i> « ~^M 2 7, (C28) 



and 

Ci Piy h k 2 
6 l + /3a 2 



( C29 ) 



The latter equation shows that the brane bending can be neglected on super-horizon scales. Eq. (|C17[) determines 
the fluctuation in the scalar field as 



.2 s 



H 



K 



Thus from Eq. (|66[) we find 



where 



—=-/*(! + /3)Ci. (C30) 



* = 47rGeff^^, (C31) 



87rG ff := n fi, 



y/1 + g 2 //i 2 
1 + /3 + 2/3H 2 /ti 2 



(C32) 



By this normalization of the gravitational constant, the relation (|C31[) gives the same result as in standard 4D 
cosmology [32]. For H 2 <C H 2 //3, we have 

G off ~ G^l + H 2 /^ 2 . (C33) 
This is nothing but the relation given in [2(| in the RS model. In the opposite limit, H 2 /i 2 / '(3, we have 

Geff * G cff . (C34) 
This agrees with what we have shown in Sec. IIVBI and thus with the result in Ref. [HI]. 



[1] D. Langlois , Prog. Th eor. Phys. Suppl. 148, 181 (2003) larXiv:hep-th/0209261| , R. Maartens, Living Rev. Rel. 7, 7 ( 2004) 
[arXiv:gr-qc/03 12059], P. Brax, C . van de Bruck and A. C. Davis, Rept. Prog. Phys. 67, 2183 (2004) [arXiv:hep-th/0404011] , 
C. Csaki, |arXiv:hep-ph/0404096| 



17 



[2] L. Randall and R. Sundrum, Phys. Rev. Lett. 83, 3370 (1999) arXiv:hep-ph/9905221 . 
[3] L. Randall and R. Sundru m, Phys. Rev. Lett. 83 , 4690 (1999) [arXiv:hep-th/9906064| . 
[4] N. Deruelle and J. Madore, [afXiv:gr-qc/0305004| 

[5] J. E. Kim, B. Kyae and H. M. Lee, Nucl. Phys. B 582, 296 (2000) [Err atum-ibid. B 591, 587 (2000)] [arXiv:hep-th/0004005l , 
J. E. Kim and H. M. Lee, Nucl. Phys. B 602, 346(2001) [arXiv:hep-th/001 0093 , K. A. Meissner a nd M. Olechowski, Phys . 
Rev. D 65, 064017 (2002) [arXiv:hep-th/0106203| , I. P. Neupane, Phys. Let t. B 512, 137 (2001) [arXiv:hep-th/0104226j , 
I. P. Neupane, Class. Qu ant. Gray. 19, 5507~[20 02) [arXiv:hep-th/0106100 , Y. M. Cho and I. P. Neupane, Int. J. Mod. 
Phys. A 18, 2703 (2003) [arXiv:hep-th/0112227j 



[6] N. Deruelle and M. Sasaki, Prog. Theor. Phys. 110, 441 (2003) arXiv:gr-qc/0306032 . 
[7] S. C. Davis, arXiv:hep th/0402151| 

[8] K. i. Maeda and T. Torii, Phys. Rev. D 69, 024002 ( 2004) [arXiv:hep-th/0309152j . 
[9] T. Kobayashi, T. Shiromizu and N. Deruelle, |arXiv:hep-th/ 0608166 
[10] T. Shiromizu, K. i. Maeda and M. Sasaki, Phys. Rev. D 62 024012 (2000) [arXiv:gr-qc/9910076j . 

[11] N. Deruelle and T. Dolezel, Phy s Rev. D62, 103502 (20 00) [arXiv:gr-qc/0004021| , S. Nojiri, S. D. Odintsov and S. Ogushi, 
Phys. Rev. D65, 02352 1(2002) [arXiv:hep-th/0108172j , C. Germani, C. F. Sopuerta, Phys Rev. Lett. 88, 231101 (2002) 
[arXiv:hep-th /0202060| , J. E. Lidsey, S. Nojiri and S. D. Odintsov, JHEP 0206, 026 (20 02) [arXiv:hep-th/0202198l , S. No- 
jiri, S. D. Odintsov and S. Ogushi, Int. J. Mod. Phys. A 17, 4809 (2002) arXiv:hep-th/02 05187| ,"PrBinetruy, C. Charmousis, 
S. C. Davis and J. F. Dufa ux, Phys. Lett. B 544, 18 3 (2002) [arXiv:hep- th/0206089], J. E. Lidsey and N. J. Nunes, Phys. 
Rev. D 67, 103510 (2003) [arXiv:astro-ph/0303168] , S. Tsujikawa, M. Sami and r' Maartens, P hys. Rev. D 70, 63525 
(2004) [arXiv:astro-ph/0406078| , K. s. Aoyanagi and K. i. Maeda, Phys. Rev. D 70, 1235 06 (2004) [arXiv:hep-th/040 8008 , 
M. Sami, N. Savchenko and A. Toporensky, Phys. Rev. D 70, 123528 (2004) [arXiv:hep-th/0408140| , T. KobayashiTOeii. 
Rel. Grav. 37, 1869 (2005) [arXiv:gr-qc/0504027| , G. L. Alberghi, D. Bombardelli, R. Casadio and A. Tronconi, Phys. 
Rev. D 72, 025005 (2005) arXiv:hep-ph/0503060 , G. L. Alberghi and A. Tronconi, Phys. Rev. D 73, 027702 (2006) 
[arXiv:hep-ph/0510267| . 

[12] C. Charmousis and J. F. Dufaux, Class. Quant. Grav. 19, 4671 (200 2) [arXiv:hep-t h/0202107 . 

[13] M. Minamitsuji and M. Sasaki, Prog. Theor. Phys. 112, 451 (2004) |arXiv:hep-th/0404166| 

[14] J. F. Dufaux, J. E. Lidsey, R. Maartens and M. S ami, Phys. Rev. D 70, 83525 (2004) |arXiv:hep-th/0404161| . 

[15] S. Muk ohyama, Phys. Rev. D 62, 084015 (2000) [arXiv:hep-th /0004067] , S. Mukohyama, Class. Quant. Grav. 17, 4777 

(200 0) [arXiv:hep-th/0006146] , S. Mukohyama, Phys. Rev. D 64, 064006 (2001) [Erratum-ibid. D 66, 049902 (2002)] 

|arXiv:hep-th7 0104185 . 

[16] H. A. Bridgman, K. a! Malik and D. Wands, Ph ys. Rev. D 65, 043502 (2002) [arXiv:astro-ph/0107245] . 

[17] C. Deffayet, Phys. Re v. D 66, 103504 (2002) | arXiv:hep-th /0205084| , C. Deffayet Phys. Rev. D 71, 23520 (2005) 
[arXiv:hep-th/0409302| , C. Gordon and R. Maartens, Phys. Rev. D 63, 044022 (20 01) [arXiv:hep-th/0009010| , B. Gumjud- 
pai, R. Maartens and C. Gordon, Class. Quant. Grav. 20, 3295 (2003) arXiv:gr-qc/0304067 , K. Koyama and J. Soda, 
Phys. Rev. D 62, 123 502 (2000) [arXiv:hep-th /00052391 , K. Koyama and J. Soda, Phys. R ev. D 65, 023514 (2002) 
[arXiv:hep-th/0108003| , K Koyama, Phys. Rev. D 66, 084003 (2002) [arXiv: gr-qc /0204047| , K. Koyama, Phys. Rev. 
Lett. 91, 221301 (200 3) [arXiv:astro-ph/0303108] , K. Koyama, S. M izuno and D. Wands, JCAP 0508, 009 (2005) 
[arXiv:hep-th/0506102| D. Langlois, Phys. R ev. D 62, 126012 (2000) arXiv:hep-th /0005025| , D. Langlois, Phys. Rev. 
Lett. 86, 2212 (2001) [arXiv: hep-th/0 010063| , D. Langlois, R. Maartens, M. Sasaki and D. Wands, Phys. Rev. D 63, 
084009 (2001) [arXiv:hep-th/0012044| , R. Maartens, D. Wands, B. A. Bassett and I. Heard, Phys. Rev. D 62, 041301 
(2000) [arXiv:he p-ph/9912464 , C. van de Brack, M. Dorca, R. H. Brandenberger and A. Lukas, Phys. Rev. D 62, 123515 
(2000) |arXiv:h ep-th/0005032 , B. Leong, P. Dunsby, A. Challinor and A. Lasenby, Phys. Rev. D 65, 104012 (2002) 
|arXiv:gr-qc/0111033| . ~ 

[18] H. Kodama, A. Ishibashi and O. Seto , Phys. Rev. D 62, 064022 (2000) [arXiv:hep -th/0004l60] . 
[19] K. Koyama, JCAP 0409, 010 (2004) [arXiv:astro-ph/0407263l . 

[20] K. Koyama, D. Langlois, R. Maartens and D. Wands, JCAP 0411, 002 (200 4) [arXiv:hep-th/0408222| . 

[21] T. Kobayashi and T. Tana ka, Phys. Rev. D 71, 12 4028 (2005) [arXiv:hep-th/0505065| , T. Kobayas hi and T. Tanaka, Phys . 

Rev. D 73, 044005 (2006) [arXiv:hep-th/0511186| , T. Kobayashi, Phys. Rev. D 73, 12 4031 (2006) [arXiv:hep-th/0602168] . 
[22] T. Hiramatsu, K. Koyama and A. Taruya, Phys. Lett. B 578, 269 (2004) [arXiv:hep-th/0308072] , T. Hiramatsu, K. Koyama 

an d A. Taruya, Phys. Lett. B 609, 133 (2005) [arXiv:hep-th/0410247] , T. Hiramatsu, Phys. Rev. D 73, 084008 (2006) 

[arXiv:hep-th/0601105| . 
[23] S. S. Seahra, Phys. Rev. D 74, 0440 10 (2006 ) [arXiv:hep-th/0602194| . 
[24] T. Hiramatsu and K. Koyama, |arXiv:hep-th /0607068 

[25] G. W. Gibbons and S. W. Hawkin g Phys. Rev. D 15, 2 752 (1977), R. C. Myers, Phys. Rev. D 36, 392 (1987), S. C. Davis, 
Phys. Rev. D 67, 024030 (2003) [arXiv:hep-th7020 8205 , E. Gravanis and S. Willison, Phys. Lett. B 562, 118 (2003) 
[arXiv:hep-th/0209076| . 

[26] P. Binetruy, C. Deffayet and D. Langlois. Nucl. Phys. B 565, 269 (2000) arXiv:hep-th/9905012 , P. Binetruy, C. Deffayet, 
U. Ellwange r and D. Langlois, Phy s. Lett. B 477, 285 (2000) [arXiv:hep-th/9910219|, S. Mu kohyama, Phys. Lett. B 473, 
241 (2000) arXiv:hep-th/9911165 , D. Ida, JHEP 0009, 014 (2000) arXiv:gr-qc/9912002 , P. Kraus, JHEP 9912, 011 
(1999) [arXiv:hep-th/9910149| . ~ 

[27] G. R. Dvali, G. Gabadadze and M. Porrati, Phys. Le tt. B 4 85, 208 (2000) [arXiv:hep-th/0005016| . 

[28] K. Koyama and S. Mizuno, JCAP 0607, 013 (2006) [arXiv:gr-qc/0606056| . 

[29] M. Sasaki, Prog. Theor. Phys. 76, 1036 (1986), V. F. Mukhanov, Sov. Phys. JETP 67 (1988) 1297 [Zh. Eksp. Teor. Fiz. 
94N7 (1988) 1], V. F. Mukhanov, H. A. Feldman and R. H. Brandenberger, Phys. Rept. 215, 203 (1992). 



18 



[30] D. Wands, K. A. Malik, D. H. Lyth and A. R. Liddle, Phys. Rev. D 62, 043527 (2000) [arXiv:astro- ph/0003278|. 

[31] H. Yoshiguchi and K. Koyama, Phys. Rev. D 71, 04 3519 (20 05) [arXiv:hep-th/0411056| , K. Koyama, A. Meimim and 

D. Wands, Phys. Rev. D 72, 064001 (2005) [arXiv:hep-th/0504201| . 
[32] V. Mukhanov, Physical foundations of cosmology (Cambridge University Press, Cambridge, 2005). 



